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Abstract 

The topological insulating phase results from inversion of the band gap due to spin-orbit coupling 
at an odd number of time-reversal symmetric points. In Bi2Se3, this inversion occurs at the T point. 
For bulk Bi2Se3, we have analyzed the effect of arbitrary strain on the T point band gap using 
Density Functional Theory. By computing the band structure both with and without spin-orbit 
interactions, we consider the effects of strain on the gap via Coulombic interaction and spin- 
orbit interaction separately. While compressive strain acts to decrease the Coulombic gap, it also 
increases the strength of the spin-orbit interaction, increasing the inverted gap. Comparison with 
Bi2Te3 supports the conclusion that effects on both Coulombic and spin-orbit interactions are 
critical to understanding the behavior of topological insulators under strain, and we propose that 
the topological insulating phase can be effectively manipulated by inducing strain through chemical 
substitution. 

PACS numbers: 03.65.Vf,71.30.+h,62.20.-x,64.70.-p, 
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I. INTRODUCTION 



The discovery of the topological insulating (TI) phase of materials has garnered intense 
interest from the condensed matter physics community and spawned numerous investigations 
to explore the nature of its origins and effects. 1-6 These materials have an insulating gap in 
the bulk, while also possessing conducting, gapless edge or surface states that are protected 
by time-reversal symmetry. The prediction and observation of this phase have been made in 
real materials, 7-13 reinforcing its status as a topic of interest and importance. While most 
investigations have focused on fundamental physical properties, the unique properties of the 
topological insulating phase suggest several practical applications, including spintronics and 
quantum computation. 14-17 However, significant progress towards technological applications 
will require deep understanding of the dependence of the fundamental physics on material 
structure and composition. In this work, we investigate the relationship between the topo- 
logical insulating phase and the elastic properties of bismuth selenide. Bi 2 Se3 is chemically 
stable, easy to synthesize, and exhibits a robust topological phase. Combined with the 
existing theoretical and experimental studies, 11,12,18-24 it has emerged as the prototypical 
topological insulator, and is a natural choice for a preliminary investigation. There has 
been some interest in the effect of mechanical strain on topological effects, 25,26 but to our 
knowledge no systematic investigation has been performed. Here we use ab initio methods 
to evaluate the elastic properties of Bi 2 Se3, and to connect these to the properties of the 
topological insulating phase. 

II. METHODOLOGY 

Density functional theory calculations were performed using the Perdew-Burke-Ernzerhof- 
type 27 generalized gradient approximation (GGA) implemented in the Quantum-Espresso 28 
code. All atoms are represented by norm-conserving, optimized, designed nonlocal 29,30 pseu- 
dopotentials generated using the OPIUM package. 31 Pseudopotentials were constructed with 
and without spin-orbit coupling. All calculations are performed with a plane-wave cutoff 
of 50 Ry on an 8 x 8 x 8 Monkhorst-Pack 32 k-point mesh. The lattice parameters were 
taken from experiments (a = 4.138 A and c = 28.64 A). 33 We calculated band structures 
along the high symmetry lines in the Brillouin-zone (BZ). After fixing the lattice parameters 
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FIG. 1. Crystal Structure of Bi2Se3 : (a) The rhombohedral crystal structure of Bi2Se3 consists 
of hexagonal planes of Bi and Se stacked on top of each other along the z-direction. A quintuple 
layer with Sel-Bi-Se2-Bi-Sel is indicated by the black square, where (1) and (2) refer to different 
lattice positions, (b) Rhombohedral unit cell of Bi2Se3. 

to their experimental values, the atomic coordinates were relaxed to generate the reference 
structure. It should be noted that the calculations give significant, nonzero stress for the 
crystal in this geometry. 

Various strains were applied relative to the reference structure, including positive and 
negative uniaxial and shear strains up to 2%, as well as several combinations thereof. The 
atomic lattice coordinates were relaxed for each strain configuration, and the total and band 
gap energies were computed both with and without spin-orbit coupling. Multiple regression 
analysis was performed to find the linear and quadratic dependence of the energy on strain 
tensor components. This yielded the elastic stiffness and stress tensors. 

It is known that in bismuth selenide the topological index distinguishing ordinary insu- 
lating from topological insulating behavior is controlled by band inversion at the T point. 
Thus, a band gap stress a T and band gap stiffness c r were defined as the linear and quadratic 
coefficients relating the V point band gap to strain, by the same procudure used to determine 
the elastic tensors. 



In both cases, only the tensor elements unique under the symmetry operations of the space 
group of bismuth selenide (R3m, No. 166) were allowed as degrees of freeedom. The stiffness 
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and stress tensors, in Voigt notation, must have the forms 
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III. RESULTS AND DISCUSSION 

The bulk crystal structure of Bi 2 Se3 is rhombohedral with space group Dg d (R3m, No. 
166) 33 , shown in Fig. 1. The primitive unit cell has two Bi and three Se atoms, and the 
atomic plane arrangement is Se(l)-Bi-Se(2)-Bi-Se(l), where Se(l) and Se(2) indicate the two 
different types of selenium atom in the crystal. In the hexagonal supercell, the structure 
can be described as quintuple layers (QL) (square region in Fig. 1 of atoms stacked along 
the trigonal axis (three- fold rotational axis). 

The band structure of Bi2Se3 and related compounds have been theoretically pre- 
dicted 12 ' 13,34,35 and experimentally observed. 10 ' 11 In our calculations, the band gap of 



Element Coefficient (GPa) 

91.8 ± 1.0 

c 33 57.4 ± 1.4 

c 44 45.8 ± 1.0 

c 66 56.2 ± 1.8 

C12 36.6 ± 1.2 

C13 38.6 ± 2.0 

cu 24.2 ± 1.8 

a[ -3.447 ± 0.007 

CTa -1.977 ± 0.010 



TABLE I. The unique elements of the elastic stiffness and stress tensors of Bi2Se3. Spin-orbit 
coupling has been included. 

the unrelaxed, experimental structure of Bi 2 Se3 is 0.3 eV, which is consistent with the 
experimental data and other calculations. 10,18,19 Fig. 2 shows the band structure with and 
without spin-orbit interaction, and both are in excellent agreement with previous results. 13 
The computed elastic and band gap tensor components are given in Table I and Table II, 
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FIG. 2. Band structure in the reference strain state of of Bi2Se3 (a) excluding (NSO) and (b) 
including (SO) spin-orbit effects. The dashed line indicates the Fermi level. 

respectively. First, we note that the gap stress shown in Table II for strain normal to the 
plane of the quintuple layers ((T3) is much higher than for strain in the plane(crf), which 
is consistent with the notion that inter-layer interactions are more important in determin- 
ing the band gap than intra-layer interactions. Second, we observe the change in sign of 
cr r when spin-orbit interactions are turned off or on. This is to be expected because the 
spin orbit interaction leads to an inversion of the conduction and valence bands. Thus, 
strain reduces the bandgap for the trivial (un-inverted) phase and increases the gap for 
the topological (inverted) phase. Third, the magnitude of the gap stress is larger when 
spin-orbit interactions are present. From a tight-binding perspective, compressive strain not 
only strengthens the Coulombic interaction between sites, increasing the associated hopping 
coefficient and reducing the conventional gap, but also magnifies the spin-orbit effect and 





NSO 


SO 


Element 


Coefficient (eV) 


Coefficient (eV) 


r r 

c n 


-67.8 ± 2.6 


35.4 ± 5.6 


c 33 


55.4 ± 3.6 


-60.1 ± 7.6 


r r 

c 44 


-58.0 ± 2.6 


23.4 ± 5.6 


c 66 


-126.6 ± 4.4 


69.4 ± 9.6 


c 12 


60.2 ± 3.2 


-33.8 ± 6.8 


c 13 


6.4 ± 4.2 


-12.6 ± 10.8 




-70.0 ± 4.2 


45.6 ± 9.0 




0.16 ± 0.017 


-1.67 ± 0.037 




4.33 ± 0.023 


-5.27 ± 0.051 



TABLE II. The unique elements of the T band gap stiffness and stress tensors excluding (NSO) or 
including (SO) spin-orbit coupling 
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FIG. 3. Band structure of Bi2Se3 near the T point as <111> uniaxial strain from 5% to 9% drives 
the topological phase transition. 



its hopping coefficient, increasing the topological gap. Thus, comparing the gap stress with 
and without spin-orbit interactions provides some insight into the effects of strain on the 
essential physics of the system. 

Using the above tensors, we can predict that the topological phase transition will occur at 
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6.4% uniaxial strain in the <111> direction. In Fig. 3, the onset of the topological insulating 
phase at 7% strain can be observed through changes in the band structure as strain increases. 
At the transition point, the Dirac cone characteristic of the phase transition is distinctly 
observable. Of course, such strains are difficult to achieve experimentally. According to the 
computed elastic tensors, around 2 GPa of uniaxial tensile stress would be required to drive 
the phase transition, well past the yield stress. However, large strains may be possible by 
introducing internal stress through chemical substitution. 
Bi2Te3 is a very similar compound to Bi2Se3, differing only in substitution of the larger 





Bi 2 Se 3 (reference) 


Bi 2 Se 3 (strained) 


Bi 2 Te3 (reference) 


Lattice parameters(A) 


a= 4.138, c = 28.64 


a = 4.358, c = 30.46 


a = 4.358, c = 30.46 


Anion radius (A) 


1.98 


1.98 


2.21 


NSO Gap (eV) 


0.02 


0.31 


0.20 


SO Gap (eV) 


0.42 


-0.06 


0.63 



TABLE III. Comparison of reference Bi 2 Se3, Bi 2 Se3 strained to match reference Bi 2 Se3 lattice, and 
reference Bi 2 Te3. Spin-Orbit(SO) and Non-Spin-Orbit(NSO) L point gaps for all three structures 
are calculated. 



tellurium in place of selenium, which increases the size of the lattice by about 6%. It is also 
a topological insulator with band inversion occurring at the L point. Given the similarity, 
one may ask if it is reasonable to view Bi 2 Te3 as intrinsically strained Bi 2 Se3. To test this 
hypothesis, we performed a comparison of Bi 2 Te 3 to strained Bi 2 Se 3 . 

In order to generate an appropriate reference structure for comparison, Bi 2 Te3 was relaxed 
under identical external stress as the reference Bi 2 Se3(this is labeled "Bi2Te3 (reference)"). 
The results are shown in Table III. The computed lattice parameters of Bi 2 Te3 are in 
good agreement with experiment. The Bi 2 Se3 lattice was then strained so that the lattice 
parameters match that of reference Bi 2 Te3(this is labeled "Bi2Se3(strained)"). Using the gap 
stiffness and stress tensors, the band gap of this strained Bi 2 Se 3 was calculated and compared 
to the computed band gap of Bi 2 Te3. Without spin-orbit interaction, the strained Bi 2 Se3 
band gap is similar to the band gap of reference Bi 2 Te3. However, with spin-orbit interaction 
turned on, the gaps are dramatically different: in the strained Bi 2 Se3, the topological gap 
closes, but in Bi 2 Te3 the topological gap is quite large. This suggests that spin-orbit effects 
are strongly dependent on the chemical identity of the anion, and that treating Bi 2 Te3 as 
strained Bi 2 Se 3 fails to capture the essential physics. 
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IV. CONCLUSION 



We have shown that strain is an important parameter for influencing the topological 
insulating phase. For bismuth selenide, the direct band gap at the gamma point, where band 
inversion occurs, responds to elastic deformation in a way that can be described by adapting 
the formalism of continuum mechanics. The critical strain at which the topological phase 
transition occurs was predicted using the derived band gap stress and stiffness tensors and 
observed in the computed band structure as a Dirac cone. While it may be possible to tune 
the band gap with external stress, more interesting is the potential for inducing strain via 
chemical substitution. Viewing bismuth telluride as chemically strained bismuth selenide, 
however, fails dramatically, hinting at a complex relationship between chemical composition, 
material structure, and the physics underlying the topological insulating phase. 
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